Vortex core dynamics induced by hole defects in antiferromagnetic nanodisks 
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Direct observation of vortex states in an antiferromagnetic layer have been recently reported [Wu, 
et al, Nature Phys. 7, 303 (2011)]. In contrast to their analogues in ferromagnetic systems, namely 
in nanomagnets, the vortex core of antiferromagnets are not expected (and have not been observed) 
to present gyrotropic or any other remarkable dynamics, even when external fields are applied. 
Using simulated annealing and spin dynamics techniques we have been able to describe a number 
of properties of such a vortex state. Besides of being in agreement with reported results, our results 
also indicate, whenever applied to antiferromagnetic nanodisks, that the presence of holes in the 
sample may induce two types of motions for this vortex. Its dynamics depends upon the relative 
separation between its core and the hole: when they are very apart the vortex core oscillates near 
the nanodisk center (its equilibrium position); while, if they are sufficiently close, the core moves 
towards the hole where it is captured and remains static. 

PACS numbers: 75.10.Hk, 75.30Ds, 75.40Gb, 75.40Mg 



I. INTRODUCTION 

Ferromagnetic (FM) disks of micron and submicron 
size are characterized by an in-plane close flux magne- 
tization, minimizing the dipolar energy. However, to 
minimizing the exchange energy at the disk center, the 
magnetization revolves against the plane, in a small 
area of only a few exchange lengths, creating the vortex 
core with two distinct polarizations, up {p = +1) or 
down {p = —1). Along with the vortex chirality, which 
accounts for the in-plane rotation sense, {q = +1) 
q = —1 for (counter-) clockwise, polarization are shown 
to govern many fundamental properties of izprtex-state 
in FM nanodisks, such as its gyrotropjipjjnodetl, switching 
of the vortex core (polarity reversal)&a and vortex-pair 
excitation. A number of technological applications has 
been proposed as long as we may control one (or both) 
of these parameters on demand. 

The scenario is expected to be quite distinct in antifer- 
romagnetic (AFM) nanodisks where the vortex appears 
as a composed pattern of two anti-aligned sub-lattices, 
presumably yielding a much more rigid structure with 
no global dynamics, as those observed for its ferromag- 
netic counterparts. Actually, indirect evidence for this 
state have beeii ^obtained by inducing FM-ordered spins 
in AFM disksLm. Only very recently, direct observa- 
tion of imprinted AFM vortex states in CoO/ Fe/Ag{001) 
disks (an AFM. layer of an AFM/FM bilayer system) has 
been reportedB. In these systems, there appears a FM- 
AFM exchange coupling, in such a way that by measuring 



certain FM spin properties it provides a probe for their 
analogues throughout the antiferromagnet, like the spin- 
flipping field, the surface order parameter, the crystalline 
anisotropy, the domain wall width etc. To our knowl- 
edge, a systematic, experimental or theoretical, analysis 
concerning vortex dynamics in AFM nanodisks is still 
lacking, although a number of works has been devoted to 
study others-aspects of such excitations in general AFM 
frame worksE3t£l. From the theoretical point of view, the 
vortex dynamic in layered antiferromagnets is very differ- 
ent whenever compared to ferromagnetic materials. For 
instance, these systems are "Lorentz invariant" and they 
bear no net magnetization. It is also expected that ap- 
plied magnetic fields causes no vortex motion in AFM 
nanodisks, in deep contrast to their ferromagnetic ana- 
logues, a fact brought about by our present results. On 
the other hand, we argue that litographically inserted 
defects may introduce a new pathway for magnetiza- 
tion dynamics, including vortex core dynamics, in AFM 
nanosy stems. This is a very interesting topic by itself, 
capable of feeding further investigation and potential ap- 
plications of these AFM systems, similarly to what has 
occured to feippmagnetic disks in presence of cavities and 
other defect SLJLJ'E£rE3. In addition, we also show that 
two types of vortex core dynamics takes place in AFM 
disks as induced by a hole. Explicitly, our results show 
that if the vortex core and the defect are quite apart, 
then the core oscillates around its equilibrium position 
(a point near to the disk center), while if they are suf- 
ficiently close then the core moves towards the cavity 
where it is captured and eventually trapped. In both 



2 



cases such a vortex core dynamics appears as a result of 
the competition between the two (AFM) sub-lattices to 
minimize exchange energy in the presence of the hole, in- 
side which exchange cost vanishes. 
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FIG. 1: Top view of an antiferromagnetic nanodisk 
with a vortex state. Namely, note the two anti-aligned 
sublattices. 




FIG. 2: Lateral view of the vortex configuration from 
Fig. ^ In the vortex core, the central region, dipoles 
develop out-of-plane projections. It must be emphasized 
that a AFM vortex consists of two subvortices^ each one 
with its own polarity and chirality. In this case, one has 
(p, q) = (+1, —1) and (—1, +1), so that as a whole, no 
net p neither q shows up. 



vanish for the whole vortex (see Figs. ^ and ^ further 
details below). This is the reason why its gyrotropic 
factor also falls off preventing dynamical manifestations, 
as the translational gyrotropic oscillatiop, iQptti^ vortex 
core observed in ferromagnetic nanodisks 



Before starting our analysis, it is useful to describe 
how a AFM- vortex can arise in nanodisks. The vortex 
state in finite magnetic samples may emerge as the 
ground-state resulting from the competition between 
exchange and magnetostatic interactions (if we focus on 
isotropic materials, anisotropic costs may be neglected). 
In these systems, such as magnetic disks, densities 
of effective magnetic charges are distributed in their 
volumes and surfaces, given by pM = —\J-M and 
(Jm = M • fi^ where n is the unit vector normal to 
the surface at every point, while M accounts for the 
magnetization. Therefore, to minimize the magneto- 
static energy the magnetization tends to point parallel 
to the dot surface. However, as one goes towards the 
center, exchange energy density increases enormously 
(in a continuum description, it is scaled with r~^, with 
r being the distance from the disk center). Then, to 
regularize exchange cost, magnetic moments tend to 
resolve against the plane, developing out-of-plane pro- 
jection, so that exactly at the center it is perpendicular 
to the disk face. Such a behavior is followed by each 
sublattice in a fifnite AFM system. Finally, it must be 
remarked that the AFM vortex state always appears as 
a composite pattern with opposite polarity and chirality 
at the sublattices, so that both parameters identically 



II. METHODS AND RESULTS 

In order to study such a state in very thin AFM disks, 
with thickness L and radius so that the aspect ratio 
L/R <^ 1, we firstly focus on each sublattice and assume 
that ^ - M ^ there (this is quite reasonable, except 
perhaps in a very small region very close to the vortex 
center). Thus, the only source of magnetostatic interac- 
tions could come from superficial magnetic charges, (Jm- 
Eventually, the vanishing of net magnetization in antifer- 
romagnets implies that magnetostatic energy also falls off 
whenever computed throughout the whole sample, pre- 
venting stability of domains in AFM systems, once typ- 
ical gain in configurational entropy is not sufficient to 
overcome domain wall energies. Now, to describe the 
system we pursue an alternative strategy substituting 
the AFM/FM bi-layer sample by a two-dimensional film 
containing an effective magnetostatic anisotropy com- 
ing ffQPii both FM and AFM exchange couplings (see 
Refs.tj'Ey for applications to FM nanodisks, even in pres- 
ence of defects). Now, the film is embeded with a regular 
square lattice defined inside a circumference of radius R 
where magnetization vector is distributed over the sites. 
In addition, the magnetostatic interactions due to the 
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presence of the magnetic charges in the lateral and top 
surfaces of the film are replaced by local potentials. The 
model is summarized by the following Hamiltonian: 

Ho = - ^Jijjli-jlj^ ^ ^K{i^k-nk,af^^h-jli, 

i,j a— 1,2 k i 

(1) 

where Jij = for sites inside the hole and Jij < for 
remaining sites of the film. Here /i^ = Mi(f)/Ms = 
ji^x + iJ(y -\- fi^z is the unit spin vector at position i {Mg 
is the saturation magnetization) and the sum z, j is per- 
formed over the nearest-neighbors spins. The term with 
positive constant, Aq,, mimics the magnetostatic energies 
at the top face of the disk, Ai, while A2 accounts for the 
lateral edge. In turn, h is the external magnetic field. 
The hole defects are introduced by removing from the 
system a number of neighbor spins around a site r^. Of 
course, the number of neighbors spins removed around a 
particular position (r^) defines the hole size Qi. The local 
unit vectors h must be normal to the surfaces at every 
point, so that fik^i = z for the disk face (in the x?/-plane), 
while hk,2 = R/\R\ = i^along the disk circumference of 
radius R. 

In our calculations we have used simulated annealing 
and spin dynamics simulations for square lattices oc- 
cupying all possible points inside a circumference with 
radius R (in most simulations we have used R = 10a, 
R = 20a, R = 25a and R = 40a; a is the lattice spacing 
constant). It is solved by employing the simulated 
annealing and the fourth-order predictor-corrector 
method. The results below have taken Ai = 0.2 and 
A2 = 2.0. Other values for Ai does not alter the essential 
physics reported here provided that < Ai < 0.28; for 
Ai > 0.28, the vortex become essentially planar and does 
not develop out-of-plane components (similar scenario 
occurs for 1.8 < A2 < 2.2). The vortex structure itself in 
an AFM nanodisk was obtained by using eq.(|^) and the 
simulated annealing techniques. The vortex configura- 
tion so emerged may be viewed as a composite structure 
of a pair of sub-vortices each of them associated with one 
sublattice. Each sub-vortex has a pattern quite similar 
to its counterpart in FM disks, but always appear with 
opposite polarity and chirality each other, yielding net 
vanishing of such quantities for the composite vortex 
(see Figs. |^ and ||). After being formed and stabilized, 
spin dynamics techniques are used to study how the 
composite vortex behaves when constant or sinusoidal, 
hext = hosm{ujt)^ external magnetic fields parallel or 
perpendicular to the plane of the disc are applied: in 
contrast to what happens to a vortex in FM disk, here 
its core remains static at the disk center (its equilibrium 
position), despite the intensity, frequency, and direction 
of the field, in agreement with theoretical expectations^ 
and experimental observation*^ . 

Dynamics has been turned on whenever defects are 
introduced, by removing some spins around a given 
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(c) The final position of vortex core (vortex trapped by the hole) 



FIG. 3: Partial view, showing only a part of the disk, 
considering the vortex core and the hole. The sequences 
of the pictures show the movement of the vortex in 
direction to the hole. 



site. Now, the hole induces translational motion of the 
vortex core yielding two very distinct consequences: 
an oscillatory motion of the vortex core near the disk 
center whenever the cavity is far away the core; or, it 
may be captured if placed sufficiently close to the hole. 
To be more precise, we have inserted a single hole of 
radius g = a (with four neighbor spins removed) in the 
system. We have chosen the origin of the coordinates 
to be the vortex center, which coincides with disk 
center, (xo^yo) = (20a, 20a) (for a lattice composed 
by 40x40 = 1600 spins), and the simulations have 
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considered the presence of a unique hole, introduced 
somewhere else. Explicitly, we have taken the hole center 
to be located at ten different points along the x-axis, say, 
at r = (a, 0), (2a, 0), . . . (10a, 0). Our results have shown 
that for a < r < 8a, the vortex core moves towards the 
hole, where it is captured yielding a practically static 
coreless vortex state (see Fig. ^). On the other hand, 
whenever the inserted hole is sufficiently apart from the 
vortex center, r >= 9a, its core develops an oscillatory 
motion around a position very close to the disk center 
(along the line joining the hole center and disk center), 
whose amplitude depends on the distance core- hole, 
diminishing as r increases. We should remark that 
intentionally inserted holes in FM disks have attracted 
a great deal of efforts in the last years and a ruij:]afe|er of 
effects on the vortex state, hap_heen predictedllSrEj and 
experimentally conffrmedll^ilSEjEj. 



III. CONCLUSIONS 

In this paper we have studied finite size circular an- 
tiferromagnetic thin films. The composite vortex state 
formed by the two AFM sublattices is show to be dy- 
namically insensitive to any applied field. However, dy- 
namics may be induced by hole defects, where vortex 
exchange energy is minimized. Once AFM vortex has 
no net gyrotropic vector its associated motion cannot 
show up (contrasting with defect induced gyrodynamics 
in FM vortex, Ref.E3); rather, its subsequent translation 
gives rise to oscillations around its equilibrium position, 
if the vortex core and defect are very apart, while if they 
are placed close, the hole captures the core, yielding to a 
static corelless vortex as the final state. Our results could 
feed further investigation on AFM nanoscaled samples 
with inserted cavities, mainly when such induced vortex 
dynamics concerns. 



ACKNOWLEDGMENTS 



The authors thank CNPq, FAPEMIG and CAPES 
(Brazilian agencies) for financial support. 



ricardo.l.silva@ufes.br 


rodrigo.dias@uf€ 


3s.br| 


apereira@ufv.bi 




winder@ufv.bi 




lucasmol@ufv. 


31 





S.-B. Choe, Y. Acremann, A. Scholl, A. Bauer, A. Doran, 

J. Stohr, and H.A. Padmore, Science 304, 420, (2004). 
2 B. Van Waeyenberge, A. Puzic, H. Stoll, K. W. Chou, T. 

Tyliszczak, R. Hentel, M. Fahnle, H. Briickl, K. Rott, G. 

Reiss, I. Neudecker, D. Weiss, C. H. Back, and G. Schiitz, 

Nature (London) 444, 461 (2006). 
^ K. Yamada, S. Kasai, Y. Nakatani, K. Kobayashi, H. 

Kohno, A. Thiaville, T. Ono, Nat. Mater. 6, 270 (2007). 
^ R. Hertel, S. Gliga, M. Fahnle, and CM. Schneider, Phys. 

Rev. Lett. 98, 117201 (2007). 
^ K.S. Buchanan, P.E. Roy, M. Grimsditch, F.Y. Fradin, 

K.Yu. Guslienko, S.D. Bader, V. Novosad, Nature Phys. 

1, 172 (2005). 

^ J. Sort, K.S. Buchanan, V. Novosad, A. Hoffmann, G. 
Salazar Alvarez, M.D. Bar, B. Dieny, and J. Nogus, Phys. 
Rev. Lett. 97, 067201 (2006). 

M. Tanase, A.K. Petford-Long, O. Heinonen, K.S. 
Buchanan. J. Sort, J. Nogus, Phys. Rev. B 79, 014436 
(2009). 

^ G. Salazar- Alvarez, J.J. Kavich, J. Sort, A. Mugarza, S. 
Stepanow, A. Potenza, H. Marchetto, S. S. Dhesi, V. Baltz, 
B. Dieny, A. Weber, L.J. Heyderman, J. Nogus, P. Gam- 
bardella, Appl. Phys. Lett. 95, 012510 (2009). 

^ J. Wu, D. Carlton, J.S. Park, Y. Meng, E. Arenholz, A. 
Doran, A.T. Young, A. Scholl, C. Hwang, H.W. Zhao, 
J.Bokor, Z.Q. Qiu , Nature Phys. 7, 303 (2011). 
^° T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, T. Ono, 
Science 289, 930 (2000). 



K.Yu. Guslienko, B.A. Ivanov, V. Novosad, Y. Otani, H. 
Shima, and K. Fukamichi, J. Appl. Phys. 91, 8037 (2002). 
^2 A.R. Volkel, F.G. Mertens, A.R. Bishop, G.M. Wysin, 
Phys. Rev. B 43, 5992 (1991). 

A. Bogdanov and A. Shestakov, Phys. Solid State, 40, 1350 
(1998). 

M.V. Chetkin, Yu. L Kurbatova, T.B. Shapaeva, JETP 
Lett. 73, 294 (2001). 

K.Y. Guslienko and A. Hoffmann, J. App. Phys. 101, 
093901 (2007). 

M. Rahm, J. Biberger, V. Umansky, and D. Weiss, J. Appl. 
Phys. 93, 7429 (2003). 

M. Rahm, R. Hollinger, V. Umansky, and D. Weiss, J. 
Appl. Phys. 95, 6708 (2004). 

M. Rahm, J. Stahl, W. Wegscheider, and D. Weiss, Appl. 
Phys. Lett. 85, 1553 (2004). 

R.L. Silva et al, Phys. Rev. B 78, 054423 (2008). 
R.L. Silva, A.R. Pereira, and W.A. Moura-Melo, J. App. 
Phys. 105, 014314 (2009). 
2^ A.R. Pereira, Phys. Rev. B 71 (2005) 224404; J. App. 
Phys. 97 (2005) 094303. 

A.R. Pereira, A.R. Moura, W.A. Moura-Melo, D.F. 
Carneiro, S.A. Leonel, and P.Z. Coura, J. App. Phys. 101 
(2007) 034310. 

W.A. Moura-Melo, A.R. Pereira, R.L. Silva, and N.M. 
Oliveira-Neto, J. App. Phys. 103 (2008) 124306. 
'^^ A.R. Pereira and W.A. Moura-Melo, ''Vortex behavior in 
ferromagnetic systems with small defects: From macro 
to nano structured magnets^\ in "Electromagnetic, Mag- 
netostatic, and Exchange-Interaction Vortices in Con- 
fined Magnetic Structures", edited by E.G. Kamenetskii, 
Transworld Research Network (2009), Kerala, India. 



5 



F.A Apolonio, W.A. Moura-Melo, F.R Crisafuli, A.R. 
Pereira, R.L. Silva, J. App. Phys. 106, 084320 (2009). 
T. Uhlig, M. Rahm, C. Dietrich, R. Hollinger, M. 
Heumann, D. Weiss, and J. Zweck, Phys. Rev. Lett. 95 
(2005) 237205. 

K. Kuepper, L. Bischoff, Ch. Akhmadahev, J. Fassben- 
der, H. Stoll, K.W. Chou, A. Puzic, K. Fauth, D. Dol- 
gos, G. Schiitz, B. Van Waeyenberge, T. Tyhszczak, I. 



Neudecker, G. Woltersdorf, and C.H. Back, App. Phys. 
Lett. 90, 062506 (2007). 

R.L. Compton and P.A. Crowell, Phys. Rev. Lett. 97, 
137202 (2005). 

F. Hoffmann, G. Woltersdorf, A.N. Slavin, V.S. Tiberke- 
vich, A. Bischof, D. Weiss, C.H. Back, Phys. Rev. B 76, 
014416 (2007). 



